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SUI-ttlARY 



A method of conformal transformation is developed 
that maps an airfoil into a straight line, the line 
being chosen es the extended chord line of the airfoil. 
The mapping Is accomplished by operating directly v/lth 
the airfoil ordlnates. Tho absence of any prellulnary 
transformation is foijnd to shorten the v.'ork substantially 
over that of previous nethods. Use is made of the 
superposition of solutions to obtain a rigorous counter- 
part of the approximate methods of thln-alrfoll theory. 
The method is applied to the solution of the direct and 
inverse problems for arbitrary airfoils and pressure 
distributions. Numerical examples are given. Appli- 
cations to more general types of regions, in particular 
to biplanes and to cascades of airfoils, are indicated. 



INTRODUCTION 



In an attempt to set up an efficient numerical method 
for finding the potential flow through an arbltj^ary cas- 
cade of airfoils (reference 1) a method of conformal 
transformation was developed that was found to apply to 
advantage In the case of isoleted airfoils. 

The method consists in transforming the Isolated 
airfoil directly to a straight line, namely, the extended 
chord line of the airfoil. The absence of the hitherto 
usual preliminary transformation of the airfoil into a 
near circle makes for a decided simplification of concept 
and procedure . 



2 



MCA ARB No. l4K22a 



a?he exposition of the method, followed by its appli- 
cation to the direct problem of the conformal mapping of 
given airfoils, Is given In part I of this paper. In 
part II the method Is applied to the inverse problem of 
airfoil theory; namely, the derivation of an airfoil sec- 
tion to satisfy a prescribed velocity distribution. A 
coTiparlson with previous Inverse methods Is made. Addi- 
tional material that will be of use In the application of 
the inothod is given in the appendixes. In appendix A cer- 
tain numerical details of the calculations are discussed. 
In appendix B extensions of the method to the conformal 
mapolng of other types of regions are Indicated. The 
relation of the methods used for the mapping of airfoils 
to the Cauchy Integral formula is discussed In appendix C. 

Acknowledgment is made to ?Trs. I.ois Evans Doran of 
the computing staff of the Langley full-scale tunnel for 
her assistance in making the calculations. 



SYMEOLS 



z = X + ly plane of airfoil 

S = S Itj plane of straight lines 

p plane of unit circle 

CP central angle of circle 

Ax component of Cartesian mapping function (CMF) 

parallel to chord 

Ay component of Cartesian mapping function perpen- 

dicular to chord 

AXq, Ay© particular CMP's, tables l and TI 

T displacement constant for locating airfoil 

r = 2r diameter of circle, senilenpth of stral^t line 

Cji = aji + Ibji coefficients of series for" CKF 

P» nepatlve of central angle of circle, corresponding 

' to leading edge of airfoil 
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Pm central angle of circle minus l80°, corresponding 

to trailing edge of airfoil - 

c airfoil chord 

c^ section lift coefficient 

velocity at surface of airfoil, fraction of free- 
stream velocity 

Vp velocity at surface of circle, fraction of free- 

stream velocity 

V free-stream velocity 

ds element of length on airfoil 

r circulation 

u^ thickness factor 

Uq camber factor 

T thickness ratio 

\ normalizing constant 

k denominator of equation (17) 

C camber, percent 

6x, 6y incremental CIi9F<s 

U positive area under approximate Vp((p) curve 

L negative area under approximate Vp((p) curve 

a angle of attack 

a J ideal angle of attack 

Y = a + P-T 

true potential 
®a approximate potential 

9 central angle of near circle 

e = (p - 8 
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Subscripts : 

N leading edge (nose) 

T trailing edge 

c camber 
t thickness 

o, 1, 2 saccessive approximation In direct or Inverse 
CVF nethods 



I - THE DIRECT POTENTIAL PROBLEM OP AIRFOIL THEORY 



TKS CARTESIAN J4APPING PUIICTION 
The Derivation of the Cartesian flapping Function 



Consider the tranaf omatlon of an airfoil, z-plane. 
Into a straight line, S-plane (fig. 1). The vector 
distance between confomally corresponding points such 
as P^ and P^ on the tv/o contours Is composed of a 
horizontal displacement Ax and a vertical displace- 
ment Ay. The quantity Ax + 1 Ay Is only another way 
of writing the analytic function z - t; that Is, 

z - £ = (x + ly) - + In) 

= (x - £) + l(y - n) 

5 Ax + 1 Ay (1) 

By Rlemann's basic existence theorem on conformal 

mapping, the function z - i connecting conformally 
corresponding points in the z- and S-planes is a regular 
function of either z or i everywhere outside the 
airfoil or straight line. This function will be referred 
to as a Cartesian mappinf^ function, or CIIF. In order to 
map an airfoil onto a straight line, the airfoil ordl- 
nates Ay are regarded as the imaginary part of an 
analytic function on the straight line and the problem 
reduces to the calculation of the real part Ax. 
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The. .calculatijon or the real part .qf an. analytic 
function on a closed contour from the known values of 
the Imaginary part Is well known. It Is convenient for 
this calculation to consider the straight line as con- 
formally related to a circle, p-plene, by the familiar 
transformation 

£ • T = p + si (2a) 
P 

where the constant displacement t has been Inserted for 
future convenience In locating the airfoil. For corre- 
sponding points on the stralgtit line and the circle, 
equation (2a) reduces to 



£ = T + r cos o 
r = 0 



(2b) 



Considered as a function of p, therefore, the CI-ll'' z - t 
la re/jular everTvvhero outside the circle and Is therefore 
expressible by the Inverse power series ? 



"S=^^^ (5) 



0 p" 



The analogy of equation (5) with the Theodorsen-Crarrlck 
transformation (reference 2) 



which relates conformally a near circle, pi -plane, to a 
circle, p-plane, may be noted. On the circle proper, 

where p = Re^*"*, and defining Cn = ^ "*' ^^nt equa- 
tion (3) reduces to two conjugate Fourier series for the 
CMF; namely. 



00 



Ax = a +XI COS n<? + 2Z si-" i' (U) 
1 1 R^ 



1 
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Ay = bo + > cos ncp - sin ncp 

1 r" 1 R^ 



(5) 



These series evidently determine Ax from Ay or vice 
versa. 

An alternative method of performing this calculation 
is possible. It is knov/n that if the real and imaginary 
parts of a function are given by conjugate Fourier series, 
as in equations (i;) and (5)* v/lth the constant terms 
zero, two Integral relations are satisfied. (See, for 
example, references 2 and 3j also, appendix C.) These 
relations are 

Ax(cp) = - ^f^" Ay(«P') cot 2J_JLJ£ d©' (6) 
2TrJo 2 

Ay((p) = -^1 Ax(qp') cot ~ ^ dtp' (?) 

Before tho detailed application of the CJliF z - i 
to the solution of the direct and inverse problems of 
airfoil theory is made, some necessary basic properties 
of this function will be discussed. 



Airfoil Position for Given CLIP 



It is noted first that tho regions at infinity in 
the throe planes are the same except for a trivial and 
arbitrary translation; that is, by equations (1), (2a>, 
and (5), 

11m z - £ s AXa, + 1 Ay„ = Cq = + ibo 
z, 

11m £ = p + T 
i» P-»« 

Secondly, if an airfoil is to be mapped into a 
straight line, it becomes necessary to know the point on 
the straight line corresponding to the trailing edge of 
the airfoil. For a given CMF, Ax(cp), Ay(«P), and 
straight line of length 2r located as in figure 1, 
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the airfoil coordinates x, y are obtained from equa- 
tions (1) and (2b) as ■ - 

X = T + r cos <? + Ax((P) (9) 

y = Ay(<p) (10) 

The leading and trailing edges of the airfoil will be 
taken as the points corresponding to the extremities of 
the airfoil abscissas. The corresponding locations on 
the circle are therefore determined by maximizing x 
with respect to (P in equation (9)> Thus 

^ = 0 = -r sin <r + 



d<p d<p 
or 

sin (P = ^ (11) 
r d<r 

The condition (11) yields (usually by graphical deter- 
mination) the angles corresponding to the leading and 
trailing edges (fig. 1) 



S ^ + ^T 



(12) 



It will be found convenient to so alter the position 
and scale of a derived airfoil that, for example, its 
chordwiae extremities are located at x = ±1 and the 
trailing edge has the ordinate y = 0 (to be referred 
to as the normal form) . The chord c of a derived air- 
foil is by definition the difference in airfoil abscissa 
extremities, or by equations (12) and (9)* 

c = r^cos P[j - cos Ptp^ + Ax^CPfyJ - Ax^tftp) (IJ) 

The increase in scale from c to some desired Cq is 
obtained simply by multiplying r. Ax, and Ay by the 
factor Cq/c. The translation necessary to bring the 
trailing edge of the airfoil to its desired location is 
then accomplished by adjusting the translation constants t 
and bg . 
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Velocity Distribution on Airfoil 

Once the CMP Ax(cp), Ay(<P ) and the diameter of" 
circle r of an airfoil have been detemlned, the 
velocity at a point on Its surface Is obtained In a well- 
known manner as the product of the known velocity at the 
corresponding point of the circle and the stretching 
factor from the circle to the alrfoU; that Is, 



vz(cp) = r §1 Vp(q)) (ll+) 

where Vp((p) Is half the velocity on the circle (since 
r = 2R) and ds Is the element of length on the airfoil, 

The velocity on the circle Vp(cp), which makes the 
point <P = TT + Pfji corresponding to the trailing edge 
of the airfoil a stagnation point (Kutta condition), is 

Vp(qp) = |sln (CP + a) + ain (a + B^^j (15) 

where a Is the angle of attack. The velocities Vp 
and Vg are expressed nondimenslonally as fractions of 
free-stream velocity. The stretching factor ds/d<p is 
obtained from equations (9) &nd (10) as 



The velocity V2(<p), equation HU.) , therefore becomes 

Isln (cp + a) + sin (a + 0m)l 
V2((P) = \ ^ (17) 

This equation is the general expression, in terns of the 
CMP, for the velocity at the surface, equations (9) and 
(10), of an arbitrary airfoil. The denominator depends 
only on the airfoil geometry, while the numerator depends 
also on the angle of attack. Equation (1?) is similar 
to the corresponding expression in the Theodorsen-Garrick 
method except for the absence of the factor representing 
a preliminary transformation from the airfoil to a near 
circle . 
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The expressions for the lift coefficient and Ideal 
Wangle of attack nay be noted. The circulation T around 
tlie airfoil is (V is free-stream velocity) 

r = RV sin (a + Pq.) (18) 

The lift coefficient c^ is defined by 

I ccjV 5 r 

Hence 

ci = Utt £ sin (a + 6rp) (19) 

where the airfoil chord c is given by equation (13)> 

The ideal angle of attack (reference 2) is defined 
as that angle of attuok for v.'hioh a stagnation point 
exists at the leading edge; that is, V2 = 0 for tp = -Bjt 
in equation (17)« Hence, * 

<ij = ^ (20) 



Superposition of Solutions 

The sun of two analytic functions is an analytic 
function; therefore, for a given p-plane circle, the sun 
of two CMF's is itself a CMF as iu also evident from 
equations (I|.) to (7). Thus, the C£!F's Lx-^ + i.Ay^^ and 

+ 1 Ay2 of two component -airfoils may, for the sane r 
be added together to tjive a CMP ^Ax^ + AX2) + ^('^yi ^^2) 
and thence, by equation (17)» an exact velocity distribu- 
tion for a resultant airfoil. The resultant profile and 
Its velocity alstrlbutlon is a superposition in this 
sense of the compjanent profiles and velocity distributions 
Thus, without sacrifice of exactness and with no great 
increase of labor, airfoils may be analyzed and synthe- 
sized in terms of component- symmetrical thlclmess distri- 
butions and mean camber lines. This result provides a 
rigorous counterpart of the well-known approximate super- 
position methods of thin-airfoil vortex and source-sink 
potential theory. 
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As a particular case of superposition, a known CMF 
Lx. + 1 Ay may be multiplied by a constant S and the 
resulting CMP S Ax + IS Ay determines a new profile by 
the new displacements S Ax, S Ay from points on the 
original straight line. It is evident that, except for 
the corrections (S - 1) Ax to the airfoil abscissas, 
this new profile is increased in thickness and camber 
over the original profile by the factor S. The effect 
on the velocity distribution is that of multiplying the 
derivatives in equfitlon (17) by S. By virtue of a reduc- 
tion in scale by the factor 1/3 this profile may also 
be regarded as obtainnd from the original one by using 
the same Ax, Ay but a length of line l/S times the 
length of the original one. 

The use of superposition as well as the application 
of the CMP to some particular airfoils will be illustrated 
next. 



Application of the CI.iP to Some Particular Airfoils 

Symmetrical thicknesa distributions .- Ttie Cartesian 
mapping fimction was calculated for a symmetrical ^0- 
percent thickness ratio Joukowski profile from the known 
conformal correspondence betv^een a Joukowski profile and 
a straight line. The Cl^IF is given in normal form in 
table I. The associated constants Tq and r^ are 

given in table II and the profile itself, as determined 
either from, the standard formulas or from equations (9) 
and (10), is shown in figure 2(a), The symmetry of the 
profile required only the calculation of Ax(<P), Ay(cp) 
for 0 ^ cp ^ l80*^. The corresponding velocity distri- 
bution (fig. 2(b)) was obtained from equation (I7) by 
use of the computed values nf the derivatives. At the 
cusped trailing edge the velocity as given by equation (17) 
is indeterminate J however, the limiting form of equa- 
tion {17)* determined by differentiation of ntamerutor 
and denominator, is 




jcos (tp + a)| 



(21) 
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It Is seen from this expression that the velocity at a 
cusped edge depends on the- second derivatives of- the 
mapping function, that Is, on the curvature at the cusp. 

The computed second derivatives d^AxQ^/d^* ^^7ot/^^ 
of the CMF of table I are plotted in figure 3 for a range 
of values of (p near l80*^. 

The CMF's for synmetrlcal profiles of different 
thickness ratios werq determined from that for the 
Joukowski profile as indicated previously in the section 
"Superposition of Solutions." The factor u*. by which 
to multiply Axg, Ayg to obtain a profile of thickness 
ratio T Is obtained from 



[ ax fa) - Ax(<Prr)] 



= T 



+ ut 



where Ayg Is the maximian airfoil ordinate of the knov/n 
CI-IF' (table I) and the denominator represents the semichord 
of the derived profile. Tlie solution for u^ is 




Values of u^ were calculated from tlixs formula for 
thickness ratios of Sij. percent and 12 percent and are 
given in table II. Tho resultin/; ClIP's were thon nor- 
malized as indicated in tho section "Airfoil Position 
for Given CMP" so that the actual factors by which to 
multiply the original AXq, Ay^ were ^u^. These values 
are given in table II, together with the associated 
constants t and r. The profiles thus determined are 
shown in figiore 2(a) and tho corresponding velocity dis- 
tributions in figure 2(b). 

Tho derived profiles' are not Joukowsld profiles. 
Tho point of maximim thickness is shifted back along the 
chord somewhat as the thickness ratio decreases. Con- 
versely, the point of naxinun thickness would be shifted 
forward by going from a thin Joukowski profile to a 
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thicker one. (This result was the reason for starting 
from a thick section.) The CMP for the 12-percent thick 
derived profile Is Illustrated In figure I4.. It Is to be 
noted that the horizontal displacement function AXq^((p) 
Is symmetrical about <P = it, whereas the vertical dis- 
placement fionctlon A7ot(9) antlsymmetrlcal about 

flP = TT. 

Mean camber line s . - The CMP was next calculated for 
a circular-arc profile of 6-percent camber from the known 
conformal correspondence between a circular arc and a 
straight line. The normalized CMP and Its derivatives 
are given In table III. The C??P is Illustrated In fig- 
ure If. The symmetry in this case is with respect to 
CP = 90° and (P = 270°, the AXqq(<P) being antisymmetrica] 

and Ayoc(tP) symmetrical. The circular-arc mean camber 
line Is shown in figure 5(a) and the corresponding 
velocity distribution in flgure"5(^)" 

Derived mean camber lines were obtained from the CMP 
for the circular arc in a manner similar to that for the 
symmetrical profiles. The expression determining the 
factor Uq for a desired percent camber C is 



max 



2 jrQ cos 
with the solution for u 



j^o cos <Pjj + Uq Ax(cp^)j 



= C 



c 



2C ro cos q^j 



u^ = ^ (25) 

^^omax " 

The angle <Pjj in equation (25) (as in equation (22)) 

corresponds to the extremity of the derived mean line. 
Because the factor Uq is to multiply the derivative 

dAXg ( cp) /dcp, the angle (P^ as determined by the maxi- 
mum condition (11) depends on Uq . One or two trials 
are sufficient to deteiTTiine Uq simultaneously with 

from equations (2$) and (11) for a given desired camber C. 
Values of u^. and ^ (also (Rp by symmetry) are 

given In table IV for derived cambers, of 3 and 9 percent. 
The actual multiplying factor to obtain the derived 
CMP's in normal form is given in table IV as ^iU_. 
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The derived camber lines are shown In figure 5(a)' 
■It Is seen that the derived- oember lines have been 
separated Into distinct upper and lower surfaces. Fur- 
thermore j for the $-Tperc&nt camber line the "lower" 
surface, that Is, the surface corresponding to the lower 
part of the straight line or circle, lies above the 
"upper" surface. Although such a camber line Is physi- 
cally meaningless by Itself, nevertheless Its CMF can- be 
compounded with that for a thickness distribution to give 
a physically real result (If the resultant profile la a 
real one). The velocity distribution of the ^-percent 
camber line Is given in figure 5(b). The "velocity dis- 
tribution" of the 9 -percent camber line is included in 
figure 5(h) for arithmetical comparison although it Is 
physically meaningless for the reason Just mentioned. 

The velocities at the cusped extremities of the 
camber lines are given by equation (21). The second 
derivatives of the CMP of table III were computed. They 

are plotted in figure 5 as d^Ax^^^/dcp^, d^AyQ^/dq)^ for 

a range of <p near l80°. These second derivatives, in 
combination with those for the symmetrical profile, can 
be used to give a more accurate determination of the 
velocity at and near a cusped trailing edge than is 
obtained by using equation (17) near the trailing edge. 

Combination of syrmetrlcal profile and mean camber 
line . - The GLlF's derived for the symmetrical profiles and 
i"or the mean camber lines can now be combined in varying 
proportions to produce airfoils having both thickness 
and camber. These airfoils may be useful in themselves 
or, as in the follov/ing sections, may be used as initial 
approximations in both the direct and Inverse processes. 

As an i^llustration of such combinations, the CMF 
of the 12-percent thick symmetrical profile of figure 2(a) 
and the CMP of the 6-percent camber circular arc of 
figure 5 (a) were added together. The airfoil profile 
thus determined is shown in figure 6(a). For comparison, 
the airfoil obtained in the manner of thin-alrfoll theory 
(see, for example, reference i^.) by superposition of the 
same symmetrical profile and a 6. 5 -percent camber cir- 
cular arc (In order to duplicate the camber of the exact 
airfoil more closely) is indicated in the figure. The 
velocity distribution of the dotted airfoil should, 
according to thin-alrfoll theory, be the sum of the 
symmetrical-profile velocity and the Increment above the 
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free-stream value of the camber-line velocity. This 
velocity distribution, determined from the two component 
exact distributions at zero angle of attack, is shown 
dotted in figure 6(b). The exact velocity distribution 
of the "exact" airfoil of figure 6(a) was determined 
for the same lift coefficient (cj = 0.88, a = 1°15») 

from the known CMP. This distribution is shown In 
figure 6(b). The two velocity distributions differ ap- 
preciably, although in the directions to be expected 
from the differences in shape of the corresponding air- 
foils. 

It appears that the CI,!F*s of a relatively small 
number of useful thickness distributions and camber lines 
would suffice to yield a large number of useful combi- 
nations of which the (perfect fluid) characteristics could 
be determined exactly and easily in the manner indicated. 

The superposition of solutions can also be used with 
the airfoil mapping methods based on the conformal trans- 
formation of a near circle to a circle. There la a 
decided advantage, hov/ever, in working with the airfoil 
ordlnates directly, both in the facility of the calcula- 
tions and in the insight that Is maintained of the rela- 
tionship betv/een an airfoil and its velocity distribution. 



THE DIRECT POTENTIAL PROBLEM FOR AIRFOILS 



The direct problem for airfoils Is that of finding 
the potential flow past a given arbitrary airfoil section 
situated In a uniform free stream. This problem can be 
solved by a CMF method of successive approximation some- 
what similar to that in reference 2. 



Method of Solution 

Suppose an airfoil to be given as in figure 6(a). 
The chord is taken as any straight line such that perpen- 
diculars drawn from Its extremities are tangent to the 
airfoil. For example, the "longest-line" chord, that is, 
the longest line that can be drawn within the airfoil, 
satisfies this definition. The x-axis is taken along 
this chord and the origin Is taken at its midpoint. 
Suppose, in addition, an initial CMF LXq and ^Jo* 



] 



NACA ARR No. lUxZ2& 



15 



.straight lime , and chordwlse translatlpn constant 

to be given such that the corresponding airfoil has the 
same chord and is similar In shape to the given alrfoll. 
(At the worst the Initial alrfoll could be the given 
chord line Itself.) 

At the chordwlse locations Xo((p) of the Initial 
alrfoll, corresponding to an evenly spaced set of up- 
values by equation (9Tf the differences 6y]^(«P) between 

the ordlnates Ay]_(cp) of the given alrfoll and Ay^Ccp) 

of the Initial alrfoll are measured. The ordinate dif- 
ferences 6y]^(cp) determine a conJ\igate set of abscissa 

corrections 5xq^((p) In accordance either with equa- 
tions (ij.) and (5) or equation (6). The details of this 
calculation are given In appendix A. 

The Initial aemllength of straight line r^ corre- 
sponding to the Initial alrfoll Is then corrected to r^^, 
and the translation constant adjusted to so 

that the use of r^^ with the first approximate CMP 
Ax-j^ = AXq + Sx^, Ay]^ = hjQ + 6j-^ yields a first approxi- 
mate alrfoll of which the chordwlse extremities coincide 
v/lth those of tho given alrfoll. This correction is 
described In detail presently. If the first approximate 
alrfoll Is not satisfactorily close to the given alrfoll, 
the procedure is repeated for a second approximate alr- 
foll, and so on. The successive airfoils thus deter- 
mined provide a very useful criterion of convergence to 
the final solution; namely, the given airfoil. Evidently, 
the fundamental relation between an alrfoll and Its 
mapping circle 




can be used in the manner indicated to effect directly 
the transformation of an alrfoll Into a circle. It 
appears preferable, howevez^ to subtract R^/p from the 
second term on the right and thence to Introduce the 

' r2 
stralght-llne variable S = p + 

The exact velocity distribution of any of the 
"approximate" airfoils (hence the approximate velocity 
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distribution of the given airfoil) may be obtained from 
equation (I?) using the derivatives of the corresponding 
CMF. The zero-lift angle Pip to be used In equation (17) 
Is debermlned for each approximate airfoil along vrlth the 
corresponding correction for r. 

Tlie correction for r is necessary because If the 
chordwlse locations of the first approximate airfoil were • 
computed by equation (9) with tlie original vr.lues of r 
and T, AX]_((p) being used instead of AXq((P), the re- 
sulting chordwlse extremities v/ould in general not be at 
X = ±1. It Is therefore necessary to adjust r^ and 
such that with the derived Ax^, Ay^, 

= 1 

where cpjj^ and cptj,^ are the angles on the circle corre- 
sponding to the extremities of the desired airfoil. This 
operation was mentioned in the section "Superposition of 
Solutions." It may be termed a horizontal stretching of 
the given airfoil. The condition given by equations . (.21).) 
applied to equation (9) yields 

1 = + ri cos cpy^ + Axi^cpjj^^ "1 

f (25) 

-1 = + ri cos «Piji^ + Axi/«Pj \ 




Subtraction of the second of these equations from the 
first gives for r^ 

A^l('Ti) - ^^i('hJ 

1 + 

2 
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Addition of eqiiatlohs (23) gives for t. 

oos (Pjjj^ + COS q)Tj^ ^ ^^(^Ni) ^^ifai) 



^1 = - 



'1 



(27) 



The angles fp^^ and (Pip^ In equations (26) and (27) 

correspond to the extremities of the desired airfoil. 
They are given by graphical solution of equation (11) 

dAxi(<p) 

sin cp = ±- — (11) 

ri dq> 

Equation (11) must bo solved simultaneously with equa- 
tion (26) for ri, (Pij^, and In practice only a 

few successive trials are neoessnry. Thence Is 

obtained by equation (27) • "^e angle (p^^ determined in 

this process is equivalent to the zero-lift angle of the 
airfoil, equation (12). 



Illustrative Example of Direct Method 

As a numerical illustration of the direct method the 
velocity distribution of the NAGA 6512 airfoil was cal- 
culated. In order to obtain an initial airfoil, the CMF 
of the 6-percent camber circular arc (tables III and IV) 
was added to the CMF of the 12 -percent thick symmetrical 
profile, derived from that of table I as Indicated In a 
previous section. Before this addition was made, the 
CMF for the circular arc was increased in scale (multi- 
plied) by 1.0928/1.0072 to correspond to the same length 
of straight line r as the symmetrical profile CMP. The 
normalized resultant CUF and the associated constants are 
given In tables V(a) and VI, respectively. The Initial 
airfoil is shown in figure 7(a). 

The given airfoil, NAGA 6512, was so rotated through 
an angle of -0.88° (nose down) as to be tangent to the 
initial airfoil at the leading edge. The convergence 
near the leading edge was thereby accelerated. The given 
airfoil is shown in this position in figure 7(&)> Two 
approximations wore then carried out In accordance with 
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the procedure given In the preceding section. The numeri- 
cal results are given in tables V and VI. The first 
approximate airfoil is indicated by the circles in fig- 
ure 7(&)s the second approximate airfoil was indistin- 
guishable to the scale used (chord = 20 in.) from the 
given airfoil. The velocity distributions of the initial, 
first, and second approximate airfoils are given In fig- 
ure 7(h), together with those corresponding to one 
approximation by the Theodorsen-Garrick method (refer- 
ence 5)» second approximation velocity distribution 
differs appreciably from that of the Theodorsen-Garrick 
method on the upper surface but agrees fairly well on the 
lower surface. The discrepancy for the rearmost 5 percent 
of chord on the lower surface appears to be due to lack 
of detail in this region in the Theodorsen-Garrick cal- 
culation. 

The convergence of the CMF method is seen to be 
rapid, considering the approximate nature of the initial 
airfoil, although two approximations are required for a 
satisfactory result. The second approximation could 
probably have been made unnecessary by suitably adjusting 
the first increment 6j^{(t>) near the leading and trailing 

edges on the upper surface before calculating 6xi(<P). 

The direction in which to adjust the increment is obtained 
by comparing the thickness of the initial airfoil with 
that of the given airfoil in these regions. Because a 
thicker section has a greater concentration of chordwlse 
locations toward the extremities, for a given set of 
Q) points, than does a thinner section, the chordwlse 
stations would be expected to be shifted outward as the 
thickness of the section is Increased. The ordlnates 
Ayi((P) should therefore have been chosen at chordwlse 
stations slightly more toward the extremities than those 
given by equation (9)> 

The accuracy of the velocities is estimated to be 
within 1 percent. It was expected, and verified by pre- 
liminary calculations, that the results would tend to 
be more Inaccurate tov/ard the extremities of the airfoil 
than near the center. This result is evident from equa- 
tion (17)- A given inaccuracy in the slopes dAx/d(p and 
dAy/d(p can produce a large error In the velocity near 
the extremities, where sin CP approaches zero. This 
disadvantage does not appear in the Theodorsen-Garrick 
method, in which sin <p is replaced by one. Excessive 
error in these regions can be avoided in various ways. 
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If the initial airfoil, for which the slopes diiX-/dcp 
and dAyo/d<P have presumably been computed accxirately, 
is a good approximation in these regions, as evidenced 
by the smallness of Ox^, 6yi compared to Ax^, Ayof 

the effect of inaccuracy of the slopes dOxj^/dcp, • dOy^/dcp 

will be reduced, since they are added to the initial 
slopes dAXg/dCp, dAy^/dxp. It was to reduce the magnitude 
of tlie Incremental CMF near the leading edge that the 
NACA 6512 airfoil was drawn tangent to the initial air- 
foil in this region. 

The error in the derivatives can also be avoided by 
computing them from the differentiated Fourier series 
for 6x-j^, dy^. (See appendix A.) This calculation was 
made in the illustrative example, after it was found that 
an error of about 5 percent in the velocity on the upper 
surface leading edge could be caused by unavoidable 
inaccuracy in measuring the incremental slopes. 

IThe fact that the computed derivatives do not repre- 
sent the derivatives of the CMP but rather the deriva- 
tive of its Fourier expansion to a finite number of 
terms may introduce Inaccuracy. (The derivative Fourier 
series converges more slowly than the original series.) 
A comparison of the computed derivatives with the measiu?ed 
slopes will Indicate the limits of error, however, as v/ell 
as the true derivative curve. 

The importance of knowing the CMF derivatives ac- 
curately may make It desirable to solve the direct 
problem from the airfoil slopes, rather than from the 
airfoil itself, as given data. This variation of 
technique enables the CMP derivatives rather than the 
CMP itself to be approximated initially. Further 
details are given in reference 1. 



II - THE INVT3RSE POTENTIAL PROBLEM OP AIRFOIL THEORY 



The Inverse potential problem of airfoil theory may 
be stated as follows: Given the velocity distribution 
as a function of percent chord or surface arc of an unknown 
airfoil - to derive the airfoil. Before the questions of 
existence and uniqueness of a solution to the problem as 
thus stated are discussed, several CMF methods of solu- 
tion will be outlined and Illustrated by numerical 
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examples. Various previous methods of solution will then 
be described briefly and their Inherent limitations and 
restrictions on the prescribed velocity distribution will 
be compared with those of the CMF methods. 

The prescribed velocity distribution is assumed to 
be either a double-valued continuous function of the 
percent chord or a single-valued continuous function of 
percent arc. (Isolated discontinuities in velocity are, 
however, at least in the percent-chord case, admissible.) 



CMF Method of Potentials 

This inverse method la based on the fact that, if 
the airfoil and its corresponding flat plate and circle 
are Immersed in the sane free-stream flows and have the 
same circulation, confomally corresponding points In 
the three planes have the sane potential. 

Consider first the case where a velocity distribu- 
tion corresponding to a symmetrical airfoil at zero lift 
is specified as a function of percent chord. If an 
Initial airfoil is assumed, the prescribed velocity can 
be Integrated along its surface to yield an approximate 
potential distribution as a fimctlon of percent chord. 
This potential Increases from zero at the leading edge 
to a maximum value at the trailing edge. Of fundamental 
Importance to the success of the method is the fact that 
this potential curve depends mainly on the prescribed 
velocity distribution and only to a much lesser extent 
on the form of the initially assumed airfoil. The chord 
line of the Initial airfoil taken as the x-axls is next 
sufficiently extended that, in the same free-stream flow 
as for the airfoil, the potential, which in this case 
Is simply V£, Increases linearly from zero at its 
leading edge~to the same maximum value at the trailing 
edge as exists for the approximate potential curve derived 
initially. Horizontal displacements Cix betv/een these 
curves are then measured as a function of the straight- 
line abscissas and, hence, as a function of the central 
angle (p of the circle corresponding to the straight 
line. These horizontal displacements &x((r), together 
with the conjugate function Ay(cp) computed therefrom 
and the length of straight line previously determined, 
constitute a CMP for an airfoil that la a first approxi- 
mation to the unknown airfoil. The approximation is 
based on the use of a more or less arbitrary initial 
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airfoil to set up the first approximate potential. The 
exact velocity distribution tff- the ddrived' first approxi- 
mate airfoil can now be computed and compared with the 
prescribed velocity. If the agreement Is not satisfac- 
torily close, the procedure is repeated, with the airfoil 
Just derived taking the place of the one initially assigned. 

The complication introduced in the general case in 
which the prescribed velocity distribution corresponds 
to an luiSTmmetrical airfoil with circulation can be 
resolved as follows: It Is convenient in this case to 
discuss the potentials in the circle plane. The pre- 
scribed velocity distribution is transferred to the circle 
plane by means of the stretching factor, presumed known, 
of tho initially assumed airfoil; that Is, equation (ll;) 
is solved for v^((P). The first approximate potential 
distribution as a function of the central an£i;le CP is 
obtained by integrating Vp((5) through a <P-range of 2Tr 

radians (around the airfoil), starting from the value 
of 5" near zero for which Vp(cp) is zero (the front 

stagnation point). This approximate potential curve has 
a riinimum value of zero at the front stagnation point, 
rises to a maximum for the value of c? near tt corre- 
sponding to the rear stagn*itlon point, then falls to a 
minimum for the final value of cp (the front stagnation 
point), which Is an angle 2Tr radians from the starting 
'O-polnt. Tlie difference between the final and the initial 
potential minimiims is a first approximation to the circu- 
lation r« 

A circle of such diameter is now derived vihich, v/ith 
this circulation and the same free-stroara flow as for the 
airfoil, yields a potential distribution (henceforth called 
true potential distribution) that has the some maximua 
and minimum values as the approximate potential curve 
Just derived. If the maximum approximate potential is 
denoted by r^U and the decrease of potential (considered 
positive) from the maximum to the final value by roL, 
where r© is the diameter of the circle corresponding 
to tho initial airfoil, the parameter y first com- 
puted from 



TT 



U 



- L 



(28) 



2(y + cot y) 



U + L 
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by means of figure 8. The desired diameter r Is then 
given by 

ro(U + L) 

r = 7-f ^ T <29) 

Z4.(cos Y + Y sin y) 

The parameter y actually the sum of the angle of 
attack and zero-lift angle of the unknown airfoil, to a 
first approximation; that Is, 

Y = a + (30) 

It Is related to the circulation T by equation (l8). 

This procedure for the calculation of the diameter 
(see, for example, reference 6) follows easily from the 
expression for the potential distribution on a circle, 
obtained by integration of equation (15) as 



*t(<P) = / Vp(cp)d<P 
N 



= rQ [cos Y"*" Y si'i Y- {JP-H- a) + (cp+ a)Bln y| (51) 



If the diameter r of the derived circle is much 
greater than the diameter rQ of the circle corresponding 
to the initial airfoil, it is desirable to Increase the 
CMP Axq, Avq of the initial airfoil by a factor suffi- 
cient to modify the initial airfoil such that it corre^ 
sponds to a circle of diameter r. A nev/ approximate and 
true potential distribution Is then obtained as described 
but by using the modified initial airfoil. 

The first approximate horizontal displacement func- 
tion Is now determined as the sura of the horizontal 
displacement Axq((P) corresponding to the (modified) 
Initial airfoil and an increment 6x2(cd) produced by 
the noncolncidence of the approximate potential distri- 
bution 4^ true potential distribubion 4^. 
This horizontal Increment may be measured between the 
two potential curves, both considered plotted against 
chordwlse position in the physical plane. With siifflclent 
accuracy this increment may be computed as the vertical 
distance between the potential curves divided by the 
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slope of the approximate .potential curve; namely, the 
prescribed velocity v^* li*, therefore, all quantities 
are considered as functions of (p 

= LXq + 5xi 

_ , <&a(<P) - ^t(^) 

" v,((p) ^52) 

The ordinate function Ay-j^(tp) conjugate to AX]^(flP) 

can now be computed and, together with Axq^CV) and the 
diameter r obtained previously, determines the first 
approximate airfoil by equations (9) and (10). Calcu- 
lation or measurement of the CMF derivatives dAX]^/d9, 
d^y]^/d(p and the use of equations (11) and (1?) then 
determine the zero lift angle and the exact velocity 

distribution of the first approximate airfoil. The angle 
of attack, to a first approximatioi^ is given by equa- 
tion (50), the value of y derived from equation (28) 
being used. This exact velocity distribution Is compared 
with that prescribed and, if the agreement is not close 
enough, the procedure can be repeated with the first 
approximate airfoil as the initial airfoil. 

In the case v/hero the prescribed velocity is speci- 
fied as a function of percent arc, then by line integra- 
tion of the prescribed velocity along the percent arc, 
the true potential distribution of the xmknown airfoil 
is known as a function of arc (except for a trivial scale 
factor). The maximum and minimum values of this potential 
distribution then permit the unique determination, by the 
calculation previously described, of the circle corre- 
sponding conformally to the unknown airfoil. Correlation 
of the potential distribution of this circle with the 
potential distribution as a function of arc initially 
calculated therefore yields exactly the potential distri- 
bution of the unknown airfoil as a function of the central 
angle q> of the circle. This fact has been noted by 
Gebelein (reference 6), The calculation of the diameter r 
as outlined above for the percent-chord case Is thus 
unnecessary. The remainder of the procedure Is the same, 
the successive approximate airfoils now being adjusted 
to correspond conformally to this circle, before corre- 
lating their percent -arc lengths with the prescribed 
velocity distribution In preparation for the next 
approxi mat ion. 
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The successive contours detepmlned by the method of 
potentials are, of necessity, closed contours, whether or 
not the sequence of contours converges to a solution 
satisfying (mathematically) the prescribed velocity dis- 
tribution. The closure of the contours Is a consequence 
of the method of setting up the horizontal displace- 
ments, &x((p), and solving for Ay((p), by which the 
contour coordinates are obtained as a ingle -valued fuiic- 
blons of 9. The necessity for closed contours does not, 
however, exclude the possibility of deriving physically 
unreal shapes; namely, contours of flguro-eight type. 
This point will be discussed at greater length later but 
It may be remarked here that it Is the extra degree of 
freedom Introduced by the class of figure-eight typo 
contours that admits the possibility of a unique solu- 
tion to the Inverse problem treated here. 

It will have been noticed that, whereas in the direct 
method a Ay is determined from the given data - that is, 
the airfoil - and a Ax is coraputed therefrom, conversely, 
in the Inverse method of potentials a Ax Is determined 
from the given data - that l3, the velocity distribution - 
and a Ay is computed therefrom. Similarly, just as 
the direct problem can also be solved by deriving dAy/dcp 
from the given airfoil slopes and thence computing 
dAx/dcp, so, conversely, can the inverse probleri be solved 
by deriving dAx/dcp fron the prescribed velocity dis- 
tribution and thence computing dAy/dCP. This Inverse 
method of derivatives will be discussed after some 
numorical examples nve presented, Illustrating the method 
of potentials. 



applie St to the derivation of the symmetrical profile 

corresponding to the prescribed velocity distribution 
shown in figure 9(a). As an initial airfoil the 12- 
percent thick profile derived from the 50-percent thick 
Joukowskl profile in part I was used. The Initial CISP 
and associated constants are given in table VII. The 
initial airfoil and its velocity distribution are shown 
In figure 9» ^li® first increment .GMF and the resultant 
first approximate airfoil and Its exact velocity distri- 
bution were calculated by the procedure of the preceding 
section. The incremental slopes d5xi/d(p, dSy^^/dcp 

were computed and found to approximate the measured slopes 



Examples of CMP Method of Potentials 




potentials was 
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very closely. The results are presented In table VII 
and, figure 5 seen that the change in velocity 

arid profile accomplished by one step' of the- inverse 
process is large; that is, the convergence is rapid. 
The high velocity of the first point on the upper surface 
(cp = 15°) is due to lack of detail in the calculation. 
(Twelve points on the upper surface were calculated.) 
For practical piuTposes the nose could be easily modified 
to reduce this velocity if desired without going through 
a complete second approximation. 

Mean camber line for uniform velocity increment . - 
As a second example of the Inverse CMP method, the profile 
producing uniform equal and opposite velocity increments 
on upper and lower surfaces was derived. By the methods 
of thin-airfoil theory this velocity distribution yields 
the so-called logarithmic camber line. The prescribed 
velocity distribution is indicated in figure 10(a). The 
velocity peaks at the extrenities of the prescribed 
velocity curve were assumed in order to compensate for 
an expected rounding off of the velocity in this region 
in working up from the initial velocity distribution. 
The convergence to the prescribed uniform velocity dis- 
tribution would thereby be accelerated. The initial 
airfoil was taken as tlie 6-percent camber circular arc, 
discussed in part I. The initial CMP and its associated 
constants are given in tables III and IV. The circular 
arc and its velocity distribution are shown in figure 10. 

A first approximation v/as calculated as outlined in 
the previous section. A numerical difficulty appeared 
In the process of solving equation (11) for the zero- 
lift angle of the first approximate airfoil. It appeared 
that a 2l|.-point calculation (12 points by symmetry) did 
not give sufficient detail In the range ir < tp < 12. it 

to yield a reliable solution of equation (11) for the 
zero -lift angle. This result was a consequence of the 
prescribed velocity discontinuity at the extremities with 
the consequent large but local changes in CMF and profile 
shape required in these regions. The solution obtained 
for the zero-lift angle was Pm = 6.1°, which by equa- 
tion (19) with r = I.OOI1.5 and a^^ = 0 yielded 
cj = 0.67- The desired c^, however, is O.8O, which 
would correspond to Prp = 7.27°. It was considered 
that a relatively minute change in the shape of the 
extremities of the derived camber line would alter the 
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slope dAxj^/dq) In the desired range sufficiently to 
yield a zero-lift angle of = 7.27°. On the other 

hand the effect of such a local change on the CMF as a 
whole would be small. The velocity distributions of the 
derived profile were therefore computed for both zero- 
lift angles quoted previously. 

The results are given in table VIII and in figure 10. 
Included for comparison in figure 10 (b) (vertical scale 
magnified) is the logarithmic mean line of thin-airfoil 
theory, computed for c^ = 0.80. The velocity distri- 
bution of the derived shape as calculated for the desired 
lift coefficient of Ct = 0.60 is seen to be a satis- 
factory approximation to the desired rectangular velocity 
distribution. The profile itself is seen to be one of 
finite thickness as compared with the single line of 
thin-airfoil theory. Airfoils Gbtained by supei*positlon 
of this type of camber line with thickness profiles would 
therefore be increased in thickness over that of the 
basic thickness form. 

The changes in velocity distribution and in shape 
of profile are again seen to be large; that is, the con- 
vergence was rapid. As is to be expected, the rapidity 
of convergence of both the direct and inverse methods in 
comparable cases is about the same. 



CMF Method of Derivatives 

Instead of approximating by the method of potentials 
to a CMP that, when differentiated, yields the prescribed 
velocity, the CMP derivatives may be obtained directly. 
The controlling equations are equations (17)f (9)» ^uid 
a modification of equation (7)* 

. Isin ((p + a) + sin fa + 
Vz((P) = ^ (17) 



dAy = A f ^ 
dp " 2Tr d<p' 



cot ^' ~ ^ d(p» 



(7a) 
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— = cos CP + •= Ax((p} 



r r 



(9) 



These equations, together with the avixlliary equations (11) 
and (l8), constitute a set of s Izcultaneous equations from 
which the CMP derivative dAx/dcp' may be detenalned from 
a prescribed velocity distribution v^. The corresponding 
airfoil is determined by integration of dAx/dcp and its 
conjugate dAy/dO). 

Consider first the case where the velocity Is speci- 
fied as a function of percent arc. As explained in the 
previous section, the constants r and y °^ final 
circle corresponding to the unknown airfoil can in this 
case be determined initially. Points of equal potential 
along the arc and circle are then found, which yield v^ 
as a function of cp. The angle of attack a in equa- 
tion (1?) is taken as some reasonable value and dAx/r dcp 
determined by successive approximation. In the first 
approxiination dAy/r dcp nay, for example, correspond to 
some known CMP- Equation (1?) is then solved for 
dAx/r dcp, for which the conjugate dAy/r dcp is calcu- 
lated next and used as a be sis for a better determination 
of dAx/r dep. The airfoil corresponding to any approxi- 
mablon Is obtained bv integration of dAx/dcp and Its 
conjugate dAy/d<p. (The method of derivatives may be 
regarded as based on the use of the function 

ip ^i-S-Z-Li, This function is regular everywhere outside 
the circle p - Re , approaches zero at infinity, and 

reduces to + 1 ^ on the circle itself. ) 

Qcp dtp 

In general the dAx/d9 as determined in any approxi- 
mation will have an average value other than zero. The 
Ax((p) obtained, say, by Integration of its Fourier 
series would therefore contain a term proportional to (P 
in addition to a Fourier series. Thus, Ax(cp) would 
not be a single-valued fvinction of (p and the resulting 
contour v/ould not close. Simply subtracting the average 
value of dAx/d(p (the constant term in its Fourier series), 
however, will close the derived contour. If the method 
converges, this average value approaches zero in the suc- 
cessive approximations , 



A preliminary over-all adjustment of an initially 
chosen CldF may be desirable. Thus, if dAx^^/dcp is 
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calculated In terms of the dAy-Q/dcp of a previous approxl- 
zaation end Is found to be larger than dAXg/dp by sozne 
factor, dAy^/dcp can be multiplied by this factor and 
the calculation of dAx^^/dtp repeated. 

Although the angle of attack may be arbitrarily set 
Initially In this calculation It should be so chosen that 
the final airfoil will coincide approximately In position 
with the Initial airfoil. After each calculation of 
dAx/d(p« the zero-lift angle firj, can be calculated, 

equation (11), which thereupon fixes a, since Y=<i+^m 
Is known. 

If the prescribed velocity distribution Is speci- 
fied as a function of percent chord, v^C^) must be 
determined In the successive approximations by use of 
equation (9)« The quantity y = a + may be deter- 
mined In each approximation as In the method of potentials 
or. In physically real cases, by equation (19) > Utia 
diameter r Is so determined that the successive airfoils 
are of a standard chord length. 

It Is evident from the structure of equation (17) 
that near the airfoil extremities where sin flp — > 0, and 
in particular at the nose of the airfoil where dAy/d(p 
is comparable to dAx/d(p in magnitude, the convergence 
by this method (and by the raethod of potentials) will be 
comparatively slow. If modifications to the airfoil only 
in the immediate nelgliborhood of the nose are required, 
it may be more expedient to apply a preliminary Joukov/ski 
transformation, that is, to use these methods with the 
Theodorsen-Garrlck transformation. 

An example of the use of the CMP method of deriva- 
tives to solve an Inverse problem is given in reference 1 
for the case of a cascade of airfoils. 



Method of Betz 

In the inverse method of Betz . (reference 7) an air- 
foil and its velocity distribution are assumed known 
(fig. 11) and a desired velocity is specified as a func- 
tion of percent arc. The new velocity and length of arc 
are specified in such a way that the extremities of 
potential are the same as on the known airfoil. ^ Both 
known and imknown airfoils then transform into the same 
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circle and. In particular, the velocities at points of 
equal potential on the two profiles can be found. 

In order to determine the profile corresponding to 
the new velocity, the complex displacement Z2 ~ ^l 
between points of equal potential on the two profiles Is 
expressed as a function of the corresponding complex 
velocities (denoted by Vg) thus. 



a \ dzg . dw/dz, 

dzT (^=2 - ^i; - - 1 - diTdiJ " ^ = vZ" 



- 1 



Hence 



Z2 



- Zi 



dzi 



(35) 



where the integration Is carried out along the known pro- 
file from the trailing edge, which is taken as coincident 
for the two elrfolls, to the point z-^. The complex 
function ^Z]^/^Z2 determined approximately from the 



known ratio 



corresponding to the points of equal 



potential by tlie argument that. Inasmuch as the two pro- 
files have nearly the same slope at corresponding points. 



the real part of 



^Z2 



- 1 la given by 



- 1. (This 



assumption, like the approximations In the GMF methods. 
Is least valid at the nose of the airfoil. The function 
Z2 - Z]^ la In fact a Cartesian mapping function.) The 

imaginary part Is then computed as the conjugate function, 
equation (7)« 

In addition to the restrictions on the velocity dis- 
tribution mentioned Initially, further conditions must 
be met In this method. If closed contours are to be ob- 
tained. Thus, the condition for closure of contour. 



i3k) 
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and the required coincidence of and at 

Infinity, lead to the following three restrictions on 
the real part R((p} of the Integrand In equation (3^) 
considered as a function of q> In the circle plane, 

J[)2tt nSrr nZ-n 

R(cp)dflp= / R(<p) cos <Pd(p= 1 R((p) sin <pd<P=0 (35) 
0 Jo Jo 



Method of Welnlg and Gebeleln 



The method of Welnlg and Gebeleln (reference 6) zaay 
be described essentially as follows: The given data are 
the same as In tlie Betz method. Consider the function 



log =~ = log 



(36) 



where 



slopes at corresponding 



points of the two airfoils (fig. 11). Since |v2;2j 
l^^ll known functions of with the data as given. 



and since log 



is regular outside the circle. 



~ ^zi ^® calculated as the function conjugate to 



log 



The 



angle being known, thereby 

determined and hence, by simple integration, the unknown 
airfoil coordinates are obtained. 

As in the Betz method, the condition for closure of 
the desired contour 



^c^^ Jc 



dw/dp 
dw/dz 



dp 



Jc 



dp = 0 



(57) 



leads to the additional restrictions on the prescribed 
velocity distribution, 
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1 0^" 

~ J -iog |.Va(a>)| d(p = ,0 ... 



1 P^'' 

- I log |v2(9)| sin (p dcp = -sin 

^ Jo 

1 

^ / log 1^2(9)1 COS <p d-p = -Tr(l - cos 2y) 
Jo ^ 

where y ^s given by equation (50). 



(58) 



Discussion of the Various Inverse Methods 



methods of Betz and of Welnlg-Gobeleln may be 
narrower in sco-po than the CKP nethods. The use 
functions such as In 



The 
somewhat 

of mapplns functions such as In equations (53) and {56) 
Is based on the ability to specify dZ2/dzi unambiguously 

in the correspondlnt» regions. This requirement appears 
to restrict the contours obtainable by these nethods to 
those boondlng slTjply connRcted regions, l-'urther investi- 
gation of this point is necessary, however. By the CIIP 
methods, f Igure-elglit contours Lave arisen In the coiu'se 
of solution of both the direct and the Inverse problems. 
(Ses the 9-percent camber derived moan line (fig. 5(a)) 
and the Illustrative examples in reference 1.) Such con- 
tours were first enoount'ared as preliminary results 
(unpublished) in using the method of potentials with the 
Theodorsen-Garrlck transformation. The CMP apparently 
maKes no fundariiental mathematical distinction between simply 
connected and figure -eight contours, for although z - £ 
must be a single -valued function of z, J, or p, the 
coordinate z itself Is of the same character as t and 
the latter has two Rlemann sheets at Its disposal In 
consequence of the Joukowskl transformation from the £- 
to the p-plane. 



The methods of Betz and of V7elnlg-0ebeleln require 
the numerically difficult closure conditions fequatlons (55) 
and (56)) to be satisfied In advance. If the nethods are 
worked through for prescribed velocity distributions 
which do not satisfy theso conditions. It appears that 
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open contoups result. In the CMP methods, . however, there 
is either no closure condition (method of potentials) or 
a numerically simple one (method of derivatives); 

r ^ ' - r ^ " = » 

[This simple closure condition in the method of deriva- 
tives is fundamentally a consequence of the fact that 
the required absence of the constant term in the inverse 
power series for the CI-!P derivative^ mapping fimction 

^ Ip 5X£_r_Ll, mentioned previously j automatically ex- 

eludes the inverse first power (the residue term) from 
the power series for d(z - O/dpH Thus, physically 
impossible velocity distributions lead to open contours 
in the Betz-Weinlg-Gebelein methods and to figure-eight 
contours in the CMP methods (if the latter converge). 
Prom the practical point of view In these cases, it may 
be easier to obtain the airfoil corresponding to the 
"best possible" physically attainable velocity distri- 
bution by tlie CMF methods than by the others. If the 
succession of airfoils deternilned by an inverse CMP method 
is seen to tend toward the development of a figure-eight, 
the successive approxiraations can be stopped at the "best 
possible" physically real airfoil. 

As to the existence and uniqueness of a solution to 
the Inverse problem as stated, a rigorous discussion of 
the solutions, for a prescribed velocity distribution, 
of the controlling equations (I7)* (7^), and (9) is 
lacking. Por physically possible velocity distributions, 
however, specified as a function of percent arc, the 
Welnig-Gebeleln method shows that there is one and only 
one airfoil as a solution. If, however, the velocity is 
specified as a function of percent chord, some further 
condition is necessary. This requirement is evident from 
the fact that one velocity distribution for an airfoil 
can, for differently chosen chords, be expressed as a 
different function of percent chord in each case. One 
chord with a given velocity as a function of percent 
chord can therefore have more than one corresponding 
airfoil. There is reason to suppose that the further 
condition for uniqueness of solution in this case is, 
the chord being defined as in the section "The Direct 
Potential Problem for Airfoils," that the ordinates to 
the airfoil at the chordwlse extremities be specified. 
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. „P?om the experience with the CMP methods gained to 
date, it la believed" thdt to a velocity " distribution 
specified as at the beginning of part II, and with the 
I'urther condition mentioned In the percent-chord case', 
there corresponds one and only one closed contour satis- 
fying the CMP system of equations . It la furthermore 
believed that the CMP methods are flexible enough to 
converge to this solution In at least those cases of 
aerodynamic interest. 



CONCLUSIONS 



■ 1. The conforraal transformation of an airfoil to 
a straight line by tlrio Cartesian mapping function (CMP) 
method results in simpler numorlcal solutions of the 
direct and inverse potential problems for airfoils than 
have been hitherto available. 

2. The use of superposition with the CMP method 
for airfoils provides a rigorous counterpart of the 
approximate methods of thin-airfoil theory. 
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APPENDIX A 



THE CALCULATION OP CONJUGATE FUNCTIONS 
BY THE RUNGE SCHEDULE 



The basic calculation for the type of mapping func- 
tion treated in this paper and in reference 2 consists 
of the computation of the real part of an analytic func- 
tion on a circle, given the imaginary part, or vice 
versa. To this end the conjiigate Fourier series, equa- 
tions ih) and (5), or the conJ\igate integral relations, 
equations (6) and (7), az*e available. This type of cal- 
culation appears to bo fundamental in many kinds of 
two-dimensional potential problems. For example, the 
solution of the Integral equation relating normal induced 
velocity to circulation in lifting-line theory can be 
solved easily by a method of successive approximation 
if the transformation from the "lifting line" to the 
circle is known. Quicker methods of calculating a func- 
tion from its conjugate than those given in this appendix 
or in reference 2 would therefore be hi^ly useful. 

The use of the Fourier series rather than the 
integral relations in the calculations of this paper was 
based on the following consideration. Because the func- 
tion l/z is regular outside the unit circle, the real 
and imaginary parts of l/z on the unit circle, namely, 
cos cp and -sin cp, satisfy the integral relations (6), 
(7)« The substitution of -sin <P for Ay in equa- 
tion (6) and subsequent numerical evaluation by the 20- 
point method of reference 2 gave results that were higher 
than cos <P by a constant error of 2.8 percent. Evalua- 
tion by a i4.0-point method reduced' the'ierror by half, or 
to l.i; percent. By the Fourier series, on the other hand, 
the first harmonic (a one-point method) suffices to give 
exact results in this case. It appears, therefore, that 
when the given real function is expressible in terms of 
a small number of harmonics, as is the case in airfoil 
applications, the Fourier series method is preferable to 
the use of the integral relations. 

The Rxmge schedule offers a convenient means of 
carrying out the basic calculation of mapping functions, 
namely, the analysis of a periodic function into its 
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Fourier aeries and the synthesis of a Poiu?ler series 
Into a function. The theory and use of the schedule Is 
desopibed, for example, in. reference .8, whe^re.l.& also 
given schedules for 12-, 2i\.~f ^6-, .and 72-point harmonic 
analyses . 

The necessaz>y analyses and syntheses in the direct 
and Inverse CMF methods are carried out In accordance 
with equations (li.) and (5) and their derivatives. 
Table IX contains the scheme of substitution into the 
Runge schedule, table X, for the various CMP methods. 
In the direct method, for example, the set of values 
6y/l2 corresponding to the evenly spaced (p-values is 
substituted into the y^, spaces at the beginning of the 
sum-table. The sums ana differences of these quantities 
are then obtained as directed at the left of the indi- 
vidual tables and substituted Into the succeeding tables. 
In this way the entire sun-table is filled out. Before 
the product-table is used, the sum-table should be checked. 

The quantities surrounded by the beavy lines in the 
svim-table are ne;:t multiplied by the proper factors at 
the left of the product-table and the results entered 
in the appropriate spaces as Indicated by the letters 
at the left of the Individual product-spaces. A heavy 
horizontal line at the lower left edge of a product- 
space Indicates that the corresponding prod\ict has 
already been obte.insd In a previous space in the sa^ne 
row. A heavy vortical line along the left edge of a 
product-space is used to enphasize that the negative 
value of the product of the sum-table quantity and the 
product-table factor is to be entered. The suma of the 
product-table columns are then entered in the I, II, III, 
and IV spaces. A check on the work of the product-table 
up to this point is provided by the columns at the rigjht. 
The sums and differences of the I, II, III, and IV quan- 
tities complete the product-table and give the Fourier 
coefficients a^, b^ corresponding to 6y. 

In order to perform a synthesis calculation from a 
set of Fourier coefficients a^, bn to the values of the 
corresponding function at the even <p-polnt3, the coef- 
ficients a^, are entered in the d and D spaces, 
respectively, in the sum- table, and the remainder of the 
siun- table and the product-table worked through as before. 
The final values in the a^, bn spaces of the product- 
table are then entered in the d and D spaces at the 
beginning of the sum-table and the sums and differences 
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obtained as Indicated by the synthesis col\mn at the 
left. (Note that do and d^^g *° ^® multiplied 

by 2.) The resulting y^^ quantities are the desired 
values of the function. 

The numerical values In tables X(a) and (b) Illus- 
trate the process of obtaining 6x]^(cp) from 5yi(<P) In 
the first approximation by the direct CMP method for the 
NACA 6512 airfoil. 
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APPENDIX B 



THE UAPPINa OF MORE GEN^L REGIONS 
Simply Cozmected Regions 



If the CMP method is applied to the mapping of a 
simply connected boundary with a vertical discontinuity, 
such as a rectangle or an infinite line with a vertical 
step, the ambiguity of the ordinate Ay at the discon- 
tinuity will prevent an automatic and rapid convergence 
of the method. Although the difficulty could be lessened 
in particular cases such as for rectangles by taking the 
diagonal as x-axis, thus romoving the vertical discon- 
tinuity, or by using symmetry, as with squares, it is 
evident that in general a reference shape particularly 
suited to the contour under investigation is needed. 
The circle has been ahown in reference 2 to be a good 
reference shape for the square . It could be expected 
therefore that an ellipse would be a good reference shape 
for the rectangle, ^''urtherraore. Just as the mapping 
function based on the circle was formed of an angular 
displacement and a radial displacement, the mapping 
function based on the ellipse should be formed of dis- 
placements along and orthogonal to the ellipse, that is, 
should be specified by elliptic coordinates. The speci- 
fication of a figure by elllpbic coordinates (\|/, 9) in 
the physical plane z is equivalent, however, to the 
transformation of the figure to a t* -plane by the tv/o 
transformations 

^ = p' ^.jr 

t» = log p' 

and specifying the transformed figure by the Cartesian 
coordinates of the t»-plane 9). The rectangle vmder 

consideration will be a near-circular shape in the p'- 
plane and a near-straight line shape in the t' -plane. 
The mapping of the rectangle by means of an elliptic 
mapping function in the physical plane is therefore seen 
to be accomplished by the Theodorsen-Garrick method in 
the near-circle p' -plane and by the CMP method in the 



\lH-ie 

where p' = e 

\ (39) 

where t' = t + 19 
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near-straight line t' -plane. Prom this point of view, 
therefore, the Theodorsen-Garrlck method consists of 
specifying an airfoil In the physical plane by elliptic 
coordinates, forming the corresponding elliptic mapping 
function - ^q) - l€, which conformal].y relates the 

airfoil to an ellipse or Joukowskl airfoil as a basic 
shape, end expressing the elliptic mapping function as a 
regular function outside the circle. On the other hand. 
In the t' = log p' -plane the Theodorsen-Garrick method 
consists of the transformation of the near-straiglit 
line ^{Q) to the straight line = Constant by means 
of what is now the ClIF - t|^o) " ^us, the 

Theodorsen-Garrick method may be regarded as a form of 
the CMF method. In which log p' takes the place of z 
and log p, the place of t. 

The mapping of simply connected regions by dif- 
ference mapping functions based on the curvilinear co- 
ordinates appropriate to the particular reference shape 
considered is therefore equivalent to using the Cf.Gi' dif- 
ference fxmction z - £ in the plane of the near-straight 
line into which the reference shape is initially trans- 
f oi>med. 



Mapping of the Entire Field 

The Fourier series representation of mapping func- 
tions, equations and (5)f enables the calculation of 
corresponding points in ths two regions to be made, once 
the correspondence of the boundaries has been calculated. 
By the latter calculation the coefficients aj^, b^ and 
the radius R of the circle of correspondence have be'en 
determined. If nov/ a larger radius R' > R be substi- 
tuted for R in equations (Ij.) and "the resulting 
synthesis of the Fourier series will yield the mapping 
function for the circle of radius R' ; that is, will 
determine points in the given plane corresponding to the 
points in the circle plane at the distance R' from the 
origin. It Is necessary, of course, to use the mapping 
function in conjunction with the shape in the physical 
plane corresponding to the larger circle. In this way the 
entire corresponding fields can be mapped out. It may 
be noted that substitution of R' < R for R in equa- 
tions (i;) and (5) enables the mapping of those corre- 
sponding points inside the original contours for which 
the resulting Fom?ler series converge . 
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It appears to bo moro difficult to find the point 
in the circle plane corresponding to a point of the given 
plane than vice versa. This calculation may, however, 
be accomplished by a method of successive approximations. 
For example, if the given plane is that of a near circle 
the polar coordinates of the given point in the near- 
circle plane are assumed to be a first approximation to 
the coordinates R' and cp of the desired point In the 
circle plane. Substitution of these values into equa- 
tions (If.) and (5) yields a first approximate mapping 
function which can be used to correct the coordinates R* 
and <D, etc. 



Biplanes 

In the case of the biplane arrangement the CMP may 
be set up directly in the physical plane in the same way 
as for the single airfoil. In place of the simple trans- 
formation from straight line to circle, however, the 
transformation from the two extended chord lines of the 
airfoils to two concentric circles is used. This trans- 
formation is derived in reference 9» The CMP method for 
biplanes bears the some relation to the nethod of ref- 
erence 9 that the CMP method for monoplane airfoils bears 
to the Ttieodorsen-Garrick method (reference 2). 

For biplanes (fig. 12) the CMP z - i, being regular 
in the region outside the two straight lines, is regular 
in the annular region of the p-plane and consequently is 
expressible as a Lavirent series in p 



- i = 



where 



2ii 

pn 



'n " ^ 



+ lb 



n 



iko) 



If, for the inner circle, the relationship Is written 



z - i = Lx.1 



1 Ayi 



p = R^e 



1<P 



ihi) 
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and for the outer circle 

z - £ = Axg + 1 

p = R2e^* 



ikz) 



there la obta5.ned, upon substitution Into equation {1^.0) 
and reduction 



Axi ( (p) = + \ ^ cos ncp + \ ^2 — ^ sin nq) (kja) 

1 1 

ro CO 

Axo(cp) = a„+\ ^'• 1 cos n9 + \ " ^"P sin n<p (l^-Jb) 

1 1 



00 



Ayi(cp)=bo + \ ''^ cos ncp -X ^" ' sin n(p (ii.3c) 



1 



CO 00 



^ b„ + b _ ^ En - a_„ 

Ay2(«P)=bo + \ -2-— a cos ncp - \ „ sin n<P (i+.5d) 

Z R2 Z_ R2 

1 1 

Those equations are the generalization to the biplane of 
equations and (5). ^e corresponding Integral rela- 
tions raay be derived as In reference 9> 

The solution of equations (1^.5) In either the direct 
or the Inverse problem may be accomplished aa before by 
succesaive approxlmatlono . x^'or example. In the direct 
method the two airfoils are given. If no initial approxi- 
mation biplane were available, the two chord lines would 
be taken as the initial straight lines. By the trans- 
formation of reference 9 this fixes the chordv/lse loca- 
tions on the straight lines corresponding to a set of 
evonly spaced cp points on the concentric circles. The 
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ordinates Ayi(cp) can therefore be measiired, which 
determines Ay2(<P) ,^7 analysis and synthesis of equa- 
tions (ij-Jc) and (i|.5d), respectively. (The radius ratio 
R2/R1 Is fixed by the Initial transformation from the 
straight lines to the concentric circles.) These Ay2((p} 
values then determine a set of Ax2(<P) values by the 
given shape of the second airfoil and the known chordwlse 
locations of Its first approximation straight line. 
Analysis of AX2((P) and subsequent synthesis of Ax-]^((P} 

by equations il^b) and (l(.^a), respectively, determines a 
correction to R^^ by a horizontal stretching process 

(constant Ax, Ay - adjustment of r-j^) to maintain the 
given airfoil chord. The procedure Is now repeated with 
^7^^^) Initial set of measured ordinates that 

determines Ay]^((P), &X]^((p), and AXp((p) as before. The 

radius R2 can now be similarly corrected. This step 
completes the first approximation. For the second approxi- 
mation a new correspondence between the corrected straight 
lines and the concentric circles Is calculated and the 
procedure repeated. 

The Inverse problem could also be solved by methods 
similar to those given for the Isolated airfoil. Sup- 
pose, for example, a wing section were given and It were 
desired to derive a slat of given chord and glvon approxi- 
mate location and having a prescribed velocity distri- 
bution. The method of surface potentials, for example, 
enables the calculation of a first approximate AX]^(tp) 

(subscript 1 refers to slat). The Initial correspondence 
of points between the straight lines and concentric 
circles, and therefore also R2/R1, being determined by 
the initially assumed straight lines, the function Ax2(<P) 

Is thereupon obtained by analysis and synthesis of equa- 
tions (I|.3a) and ([{.^b), respectively. The horizontal dis- 
placement Ax2(9) thence determines Ay2(5P) by the 
known shape of the main wing section. Tlie determination 
of Ay3^((p) by analysis and synthesis of equations (U^d) 
and (I4.5C) completes the calculation of the first approxi- 
mate slat section, for which the exact velocity distri- 
bution can now also be calculated. If the main wing 
section were also unknown then the wing section above Is 
regarded as an initial approximation, the role of the two 
airfoils Is reversed, and the procedure repeated to com- 
plete the first approximation. 
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The CtSF method can be generalized In the same manner 
for multiply connectod regions. The transformation from 
the n reference shapes (such as straight lines) to n 
circles being presumed knovm, the CMF can be set up as a 
series convergent In the region between the n circles, 
and the mapping function for each boundary explicitly 
expressed by allowing the coordinate vector to assvime Its 
value on each boundary In turn. A method of successive 
approximation for the solution of the resulting equations 
depending on the particular problem imder consideration 
would then be established. 



Cascade of Airfoils 

A simplified but practically Important n-body probleo^ 
namely, the cascade of airfoils, may be mentioned finally. 

The reference shape into v;hlch the cascade of air- 
foils, figure 15, la to bo transformed Is chosen as the 
cascade of strlght lines coinciding with the extended 
chord lines of the airfoils of the cascade. The trans- 
formation from the cascade of straight lines to a single 
circle is well-known, reference 10. The CI.!P chosen as 
indicated in figure Ip is therefore expressible as an 
Inverse power series in the circle plane and the resulting 
procedure In either the direct or the inverse problem is 
seen to be essentially the same as for Isolated airfoils. 
The detailed application of the CldF to cascades of air- 
foils is given in reference 1. 
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■ APPENDIX C 



TBE DETERMINATION OF MAPPING FDNCTIONS BY THE 



CAUCHY INTEGRAL FORMULA 



The foregoing methods of confoznnal transformation 
have been presented from the point of view of represen- 
tation of the various mapping functions as infinite 
series. In particular, the expression of the Cartesian 
mapping function as an Inverse power series valid every- 
where outside and on a circle led to the Foiirler series 
representation for the CMF on the circle itself. The 
Integral formula representation was then obtained from 
the Fourier series by the uethod of reference 5« It Is 
of interest to see how the Integral relations (6) and (7) 
can be derived directly from the Cauchy Integral formula 
for a function regular outside a circle. (These integral 
relations have also been derived by Betz, reference 7* 
by a hydrodynamlcal argunent.) Since the applicability 
of the Cauchy inte.«5ral foraiula is not restricted to 
circular boiindurles, however, trie resulba v/111 be capable 
of generalization, in principle at least, to arbitrary 
simply and multiply connected regions. 

The Cauchy integral formula gives the values of an 
analytic function f (p) within a simply connected do- 
main D in terns of its values f(t) on the boundary 
of the domain as 



where the path of Integration is counterclockwise around 
the bo\indary. Consider the domain D outside the simple 
closed boundary C in the p-plane (fig.ll].). Tills domain 
can be made simply connected by an outer boundary B and 
the cuts between the two boundaries, as indicated by the 
dotted lines. The Cauchy Integral formula for the func- 
tion f(p) at an interior point p of the domain D, 
in terms of its values on the boiondary is 




f (p) = -i- 




dt 
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whepe the equal and opposite Integrals along ^e cuts 
have been omitted. The paths of Integration are Indicated 
by the arrows In figure li|.. Tlie function f (p) Is as- 
sumed to be regular everywhere outside the boundary C 
and In particular to approach the limiting value f^, as 
p — >■<». If the boundary B Is enlarged Indefinitely, 
the Integrand of the second Integral of equation (1^.5) 
approaches f„/t and thus 

ii„ j:^ r jm dt = f„ m 



t-^-SirlJg t -p 



The variable p will now be made to approach a point t' 
on the boundary C, and equation (j+5) will consequently 
reduce to an Integral equation for the boundary values 
of a fxmctlon regular everywhere outside and on the boun- 
dary. In order to evaluate properly the contribution of 
the remaining (first) integral of equation (^5) in the 
neighborhood of t', the boundary C la modified as 
Indicated In figure li;. The point p la nade the center 
of a semicircle v^hose ends are faired into the original 
boundary. As p — >t', the modified boundory approaches 
coincidence with the original boundary. The Integral 
over the modified boimdary la now evaluated as the sum 
of liie integral over the semicircle, which in the limit. 

is half the residue of the intecrand or ^f(t')f and the 

integral over the rest of the path, which in the limit 
la analogous to the Cauchy principal value of a real 
definite Integral of which the Integrand becomes infinite 
at some point in the interval of integration. Equa- 
tion (i|.5) therefore becomes, in the limit, 

•£1*11 = J:_ f + f„ (1,7) 

2 2Trl t - t' " 

In addition, there is the auxiliary condition that 




which follows from the fact that f (p) is regular every- 
where outside the boundary C. Equation (ij.?) Is well 
known in the theory of functions of a complex variable. 
(See reference 11.) 
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If, now, the function f (p) la taken as the Cartesian 
napping f iinctlon z - £ ' or, on the boundary, 

f (t) = Ax + 1 Ay (I1.9) 

and If, further, the boundary C Is taken as a circle 
with origin at the center. 



(50) 



substitution of equations (i^9) Bl^^ (50) into equation (I|.7) 
and using equation (i|.8) (with f = 0) loads to the inte- 
gral relations (6) and (7)« I^" tii® polar mapping func- 
tion log ^ = - 16 = - - i(q) ^ 9) (reference 2) 

is substituted for f(t), the Theodorsen-Garrick integral 
relations are obtained. 



t = e^o"-"-^ 
t. = e^^"-^'' 



The Cauchy Integral formula has already been applied 
(reference 12) to problens of conformal mapping In the 
manner just indicated. Bergman has Included in refer- 
ence 12 (chapter XI) contributions of two Russian authors, 
Gershgorln and Kr^'lov. In reference 12 the mapping func- 
tion from a circle to a near circle was taken In the form 
log p. The resulting Integral equation does not appear 
to be as convenient as those of the CTilF methods . The use 

of forms such as log ^ or z - £ are not only accurate 

numerically since they express changes in the coordinates 
of the boundaries, but also they lead to pairs of integral 
equations which contain the solutions of both the direct 
and the inverse potential problems. 

Prom the analysis given it appears possible to trans- 
form conformally from one boundary to another without 
requiring the transformation from either boimdary to a 
circle, since the boundary C in equation (U?) can be 
rather arbitrary and f(t) can be taken as a mapping 
function from this boundary to another arbitrary one. 
The resulting Integral equation for the mapping function 
is, however, not as easy to solve numerically as when the 
boundary C is a circle. 
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Once the confonnal correspondence between two boun- 
daries Is known, corresponding points outside the boun- 
daries can be determined by the Cauchy integral formula 
Ikk.) ' It is noted that the Cauchy integral gives the 
correspondence of individual pairs of points rather than 
the correspondence of entire boundaries at once, which 
is given by the Fourier series representation. Further- 
more, the possibility exists of determining pairs of 
corresponding points Inside the given boundaries by the 
Cauchy Integral, that is, of analytically continuing the 
conformal transformation beyond the original domains. 
For if the transformation - from a boimdary C in a 
p-plane to a boxmdary C in a p' -plane were known, the 
outside regions corresponding, then the correspondence 
between a boundary A internal to G and a boundary A' 
internal to C, if it existed, could be determined by 
an application of the Cauchy integral formula to the 
region bounded by A and C. 

For example, if the boundaries A and C are taken 
as concentric circles and the mapping function as 



in the notation of figure 15, the Cauchy integral formula 
applied to the annular region in the p-plane (assumed 
free of singularities of the mapping function) yields, 
in the limit as the variable point p approaches the 
inner circle A, 



f (p) = log ^ 

= \Jf - i€ = (\!/ - X) - i((p - 0) 



(51) 




2"Jo ooBh (Xo -Xi) - 008 (cPo - <Pi') 



(52a) 



€i(<Pi') 



'^i *i(<Pi) cot IL-_J:- dcpi 




(52b) 
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In adt^ition, the condition of regularity of the function 
f(p) In the annular region yields the auxiliary condi- 
tions 



1 A2Tr 1 A2Tr 



1 P2Tr T /»2Tr 



O ^O 



(53) 



In the problem under consideration, the mapping function 



^o 



">^1 

e 



for the outer boundaries Is Imown. The radii » , 
of the concentric circles are given. The second Integrals 
of equations (32) are thiis knovn functions of . Equa- 
tions (52a) and (52b) therefore constitute a pair of 
integral equations, similar to those of Theodorsen-Qarrlcl^ 
for the mapping function ^if<P£) - ^^iC^l)* pertaining 
to the inner boundaries. ^ ' -^-v / 

It is noted that if the variable point p of the 
Cauchy Integral formula for the annular region is made 
to ar^proach the outer boundary C, then two additional 
Integral equations similar to equations (52a) and (52b) 
are obtained. These eqtiatlons, together with equa- 
tions (53 )» are a generalization to the case of ring 
regions of the corresponding Theodorsen-Garrick 
equations for simply connected regions and can be used 
for the oonformal mapping of near circular ring regions. 
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TABLE I 

CARTESIAN MAPPINa FDNOTION FOR S-XMHETRIGAL 
30-PERCENT lEICKNESS JODKOWSRI PROIJTLE 



i 


AXq 






dA7o/d(p' 


° ^ 12 


-0.519 


0 


0 


0.375 


1 


--30m. 


.090^. 


.119 


9 1" 

.355 


2 


Off Q 

-.250 


.lo2 


•226 


.296 


0 


- .190 


.250 


.309 


^ 

.206 


1, 


- .101 


■2o7 




.Oo9m- 


c 
P 




• 295 


.352 


-.0351 


6 


.0798 


.270 




-.lk9 


7 


.12|.8 


.218 


.212 


-.238 


8 


.169 


.150 


.0916 


-.272 


9 


.197 


.0810 


-.0261 


-.2l].0 


10 


.179 


.0291 


-.0958 


-.114.9 


11 


.I5i|. 


.OOij.12 


-.O82U 


-.0314.6 


12 


.1^2 


0 


0 


0 



TABLE II 



CONSTANTS USED WITH CMP OP TABLE I 



Profile 


T ■ 




T 


r 




(deg) 




JoTikovskl 


0.30 


1.000 


0.0887 


1.230 


0 


180 


1.000 


Derived 


.2l|. 


.805 


.0716 


1.185 


0 


180 


.835 


Derived 


.12 


.1|.02 


.0357 


1.0928 


0 


180 


.1|.53 
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TABLE III 

CMP FOR 6-PERCENT-O AMBER C IRCT3LAR-ARC PROFILE 





Slans) 






O' 


dAy^/ap 


6 


X — — 

12 


0 


0.120 


0.108 


0 


7 




.0270 


.iiij. ■ 


.0960 


-.Olj^l^. 


8 




.Oit.82 


.0953 


.0638 


-.0871 


9 




.0592 


.069I4. 


.0171 


-.109 


10 




.0565 


.oij.05 


-.0363 


-.106 


11 




.OI1.O8 


.0160 


-.08l|i4. 


-.0781 


12 




.011(2 


.00169 


-.115 


-.0279 


13 




- . 0170 


.00214.6 


-.117 


.034-6 


11^ 




-.Oil.39 


.019!]. 


-.0852 


.0926 


15 




-.0587 


.014.90 


-.0239 


.128 


16 




-.0552 


.0328 


.0506 


.123 


17 




-.0355 


.110 


.113 


.0756 


18 


! 


0 


.120 


.136 


0 



TABLE IV 

CONSTANTS USED WITH CMP OP TABLE III 



Profile 


C 

(per- 
cent) 


^c 




T 


r 


*pn 

(deg) 


(deg) 


^l 
(deg) 


°l Ideal 


Derived 


3 


0.502 


0.501 


0 


1.0052 


-3.37 


183.37 


0 


0.37 


Clrovilar 
arc 


6 


1.000 


1.000 


0 


1.0072 


-6.814- 


186.814. 


0 


.75 


Derived 


9 


1.502 


1.^+99 


0 


1.0050 


-10.26 


190.26 


0 
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TABLE V 

THE DIRECT CMP METHOD FOR THE NACA 6512 AIRFOIL* 



wATioNAL mm 

Initial Approximation™""'"^ ff'P moUmCS 



i ? X 

(radians) 

■ — 






dAJCo 








1_ 










*i 


^0 


(ci - 1.5) 




-0.145 


0.0018 


-0.126 


0.183 




0.992 


0.201 


1.633 


1 


-.168 


.0565 


-.0438 


.229 




.931 


.364 


, 1.600 


ti 


-.166 


.118 




• 236 




• o20 


«506 


1«596 


3 


-.142 


.177 


.144 


.203 




tf ptn 

.672 


tf rte 

.605 


^ tf 

1.651 


4 




.221 


.211 


.131 




• 493 


tf o e 

.685 


t sec 


5 


-.0325 


.244 


.248 


.0388 




.288 


' .742 


^ tf o c 

1.625 


6 


.0324 


.241 


.241 


-.0606 




•0678 


.783 


1.546 


rr 

7 


• Uovo 




.191 


-.100 




"•160 


.804 


1 ilOQ 

X.4iso 


Q 

O 




. Xoo 


.lU / 


— ./iLKS 




— .OOD 


.TvX 


X.<S9X 


V 


1 ^ K 




. JUOX 


-.<sXO 




"•oyy 


» Il6l 


x.xoo 


xu 




•uoox 


A*7QA 


— 




••• 




X.UOv 


11. 


.1U7 




IOC 


-.obyo 




••.9^1 


• 382 


.908 


l<e 




• UUXo 


— .XfiO 


-•UoUc 




-•993 


• 118 


DOT 


13 


•0439 


.0010 


-.0948 


.0239 




-.984 


•174 


.853 


14 


.0246 


.0094 


-.0543 


.0410 




-.894 


.452 


.834 


15 


•0156 


.0207 


-.0155 


.0452 




-.728 


.694 


.814 


16 


.0162 


.0297 


.0182 


.0245 




-.499 


.883 


.793 


17 


.0236 


.0319 


.0368 


-.0138 




-.226 


.999 


.773 


18 


.0324 


.0218 


,0252 


-.0606 




.0678 


1.024 


.759 


19 


.0337 


•0008 


-.0199 


-.0970 




.354 


• 952 


.750 


20 


•0195 


-.0263 


-.0872 


-.0987 




.555 


.792 


.742 


21 


-.0128 


-•0479 


-.151 


-.0566 




•801 


.572 


.711 


22 


-.0566 


-.0528 


-.185 


.0186 




.933 


.332 


.551 


23 


-.105 


-.0364 


-.177 


.106 




.994 


.138 


.493 



CMP' 8 of table V are with reference to chord rotated O.BSO counterclookwlM froo "long«at-line* 
chord. 



TABLE V 

THB DIRECT CKP METHOD FOR THE NACA 6512 AIRFOIL* - Contlmwd 

NATIONAL ADWSORy 



I radlana) 


1 

An 


6ji 


dxi 


A6xi 


(HP 


tx.1 


dV 




X2 


1^1 


1 

(ci = 1.5) 


0 x^S 


0,0018 


0 


-0.0057 


0.0446 


0.0151 


-0,151 


-0.0812 


0.198 


0.995 


0.192 


1.680 


1 


.055 


-.0016 


.0069 


•0429 


-.0281 


-.161 


-.0009 


.201 


.947 


.316 


1.821 


2 


.106 


-.0123 


.0138 


.0092 


-.0490 


-.153 


.0612 


.187 


.844 


.476 


1.692 


3 


.153 


-.0243 


.0113 


-.0210 


-.0383 


-.130 


.123 


.164 


.689 


.615 


1.617 


4 


.189 


-.0319 


.0024 


-.0428 


-.0203 


-.091 


.169 


.111 


.498 


.722 


1.565 


5 


.209 


-.0347 


-.0103 


-.0542 


0 


-.043 


.194 


.0388 


.278 




1 517 

X . VX 1 


6 


.213 


-.0276 


-.0248 


-.0398 


.0486 


.0076 


.201 


-.0120 


.0406 


.819 


1.475 


7 


.200 


-.0126 


-.0312 


-.0088 


.0593 


.0586 


.182 


-.0902 


-.196 


.806 


1.422 


8 


.170 


.0030 


-.0292 


.0268 


.0573 


.100 


.134 


-.150 


-.423 


.758 


1.348 


9 


.125 


.0163 


-.0179 


.0550 


.0321 


.127 


.063 


-.184 


-.627 


.672 


1.253 


10 


• 075 


.0189 


-.0068 


.0668 


-.0115 


.134 


-•0118 


-.133 


-.797 


• 524 


1.180 


11 


.030 


.0108 


.0113 


.0270 


-.0410 


.118 


-.0988 


-.131 


-.924 


.367 


1.002 


12 


.003 


.0012 


.0134 


-.0062 


-.0258 


.0864 


-.132 


-.0560 


-.993 


.129 


.828 


13 


0 


-.002 


.0120 


.0014 


-.0039 


.0559 


-.0934 


.0200 


-.986 


.176 


.835 


14 


.0057 


-.0037 


.0137 


.0076 


-.0061 


,0383 


-.0467 


.0349 


-.892 


.459 


.819 


15 


.015 


-.0057 


.0156 


.0099 


-.0088 


.0312 


-.0056 


.0344 


-.722 


.703 


.804 


16 


,020 


-.0097 


.0199 


.0145 


-.0245 


.0361 


.0327 


0 


-.487 


.895 


.782 


17 


.014 


-.0179 


.0218 


0 


-.0367 


.0454 


.0368 


-.0505 


-.210 


1.0 


.774 


18 


-.007 


-.0288 


.0192 


-.0214 


-.0429 


.0516 


.0038 


-.104 


.0846 


1.008 


.774 


19 


-.035 


-.0358 


.0067 


-.0436 


-.0065 


.0404 


-.0635 


-.104 


.361 


.914 


.786 


20 


-.059 


-.0327 


-.0030 


-.0344 


.0153 


.0165 


-.122 


-.0834 


.606 


.760 


.779 


21 


-.078 


-.0301 


-.0117 


-.0317 


.0191 


-.0245 


-.183 


-.0375 


.795 


.544 


.758 


22 


-.071 


-.0182 


-.0189 


-.0019 


-.0496 


-.0755 


-.187 


.0682 


.921 


.337 


.561 


23 


-.043 


-.0066 


-.0146 


.0262 


.0322 


-.119 


-.151 


.138 


.989 


.175 


.350 



> 
> 
> 

tz: 
o 



ro 

CO 



*CMP»a of table V are with reference to chord rotated 0,88° counterclockwise from tbngest-llne* 
chord. 



TABLE V 

THE DIRECT CMP METHOD FOR THE NACA 6512 AIRFOIL* - Concluded 



NATIONAL ADVISORY 
Second approxlmatlonCOMM'TTEE FOR AERONAUTICI 



CP 

(radians) 




072 


0X2 


dtp 


d6y2 




dAx2 

d(p 


a^j2 


^3 


^2 


V2 

VI/ J — X . o ; 


d(p 


dcp 


TI 

0 X ^ 


-0.0057 


-0.0075 


0.0043 


-0.0064 


-0.0138 


-0.146 


-0.0876 


0.184 


0.996 


0.183 


1.780 


1 


.045 


-.010 


.0005 


-.0198 


.0018 


-.160 


-.0207 


.203 


.943 


.332 


1 1.744 


2 


.098 


-.008 


-.0041 


-.0095 


.0138 


-.157 


.0517 


.201 


.836 


.488 


1.656 


3 


.149 


-.004 


-.0050 


0 


.0111 


-.135 


.123 


.175 


.681 


.617 


1.615 


4 


.188 


-.001 


-.0039 


.0092 


.0103 


-.0950 


.178 


.121 


.492 


.714 


1.584 


5 


.210 


.001 


-.0004 


.0187 


0 


-.0436 




.OOoo 


• <s7o 


#770 




6 


.212 


-.0015 


.0031 


.0069 


-.0130 


.0107 


.208 


-.0250 


.0426 


.813 


1.487 


7 


.195 


-.005 


.0032 


-.0057 


-.0130 


.0618 


.177 


-.1032 


-.193 


.812 


1.411 


8 ■ 


.162 


-.008 


.0006 


-.0160 


-.0057 


.1006 


.118 


-.1559 


-.423 


.773 


1.320 


9 


.118 


-.0075 


-.0037 


-.0141 


.0069 


.1229 


.0490 


-.1769 


-.629 


.682 


1.231 


10 


.071 


-.004 


-.0062 


0 


.0160 


.1273 


-.0118 


-.1174 


-.802 


coo 




11 


.030 


0 


-.0046 


.0133 


.0047 


.1136 


-.0855 


-.1258 


-.926 


.354 


1.030 


12 


.003 


0 


-.0020 


-.0036 


-.0033 


.0844 


-.136 


-.0593 


-.994 


.134 


.778 


13 


0 


0 


-.0015 


.0031 


-.0014 


.0544 


-.0903 


.0186 


-.985 


.178 


.840 


14 


,006 


0 


-.0006 


.0031 


0 


.0378 


-.0437 


.0349 


-.891 


.462 


.819 


15 


.015 


0 


-.0002 


.0019 


0 


.0310 


-.0037 


.0344 


-.721 


.705 


.804 


16 


.020 


0 


.0009 


.0025 


-.0027 


.0370 


.0352 


-.0027 


-.486 


.898 


.781 


17 


.013 


-.001 


.0008 


-.0008 


0 


.0462 


.0360 


-.0505 


-.209 


.999 


.774 


18 


-.009 


.001 


.0010 


.0067 


.0076 


.0526 


.0105 


-.0959 


.0845 


1.014 


.768 


19 


-.036 


-.001 


.0054 


0 


-.0191 


.0458 


-.0635 


-.123 


.365 


.915 


.782 


20 


-.064 


-.005 


.0015 


-.0076 


0 


.0180 


-.129 


-.0834 


.605 


.753 


.783 


21 


-.078 


0 


.0008 


,0015 


0 


-.0237 


-.182 


-.0375 


.792 


.545 


.751 


22 


-.073 


-.002 


.0055 


.0088 


-.0107 


-.0700 


-.178 


.0575 


.923 


.343 


.540 


23 


-.048 


-.005 


.0047 


-.0024 


-.0022 


-.115 


-.153 


.136 


.989 


.172 


.379 



*CMP»8 of table V are with reference to chord rotated 0.88° counterclockwise from ^longest-line" chord. 



TABLE VI 



CONSTANTS USED WITH CMF's OP TABLE V 



CaII angles are given with reference 'to "longest-line?' chord 
of airfoil, "Longest- line" chord is rotated clockwise 53 » 
with respect to "x-axis" chord.] 



Appr ox ima t i on 


r 


T 


Pt 




(c^ = 1.5) 


Initial 


1.1013 


O.O55I1 


6° 58' . 


7° 251 


50 33, 


1st 


1.1128 


.0550 


70 1^, 


k° 50' 




2nd 


1.1102 


.0519 


6° 55- 




50 30' 


The odor sen- 
Garrick, 
1st approximation 


1,119 




70 


70 39, 


50 21 1 
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TABLE VII.- INVERSE CMP METHOD FOR STIIMETRICAL PROFILE VELOCITY DISTRIBUTION 



> 
> 
> 

z 
o 



r- 



Initial CMP 



Plrat Increment CMF 



Flrat approximation 



9 

(radians) 



Ax„ 



ATo 



6x, 



dCXj 



AX, 



d&Xj 



d9 



0 
1 



9 

10 

11 

12 



■0.129 

-.122 
-.lOL 

-.07f3 
-.o4o5 

-.0029 
.0321 
.0597 
.0759 
.0791 

.0720 
.0619 
.0572 



.0388 

.073U 
.101 

.116 

.119 

.109 

.0877 

.0602 

.0326 

.0117 

.0017 



.oltSo 
.0909 
.12U 
.iki 

.llt2 
.122 
.085 
.036 
-.0105 
-.0585 
-.0531 
0 



0.151 

.0830 

.0559 
-.OllA 
-.0601 

-.0957 
-.110 

-.0963 
-.0597 
-.0139 

0 



0 

.0006 

-.0059 

-.oi§3 

-.02lv9 
-.0275 
-.021U 
-.0093 

-.00142 
-.0015 

-.0011 

-,0005 

0 



■..0059 
-.0156 
-.015? 
-.0071+ 
.0039 
.0151 

.Ol6lt. 
.Ollli 

.0079 
.00^6 
.0028 

0 



-.0051^ 
-.038U 
-.0389 

.OOli.3 
.01+57 

.0562 

.0099 
.0071 
.0006 

.OO^ii 
0 



.0.0151^ 

-.0539 

-.0165 
.0119 

.0370 
.01^66 
.0309 
-.0179 

-.0137 
-.0156 
-.0076 
..0091 
-.0110 



.0.1192 

-.1123 

-,1003 

-.0883 
-.0561 

-.0211 
.0200 
.0597 
.0810 
.0869 
.0802 
.0707 
.0665 



.032 

.0598 

.0867 

.1083 

.1221+ 

.1236 

.1041 

,0716 

,0405 

.0163 

.00^5 

0 



.01(26 
.0S25 

.0853 

.1^59 
.1661 
.1216 
.0W7 
.0034 
.0379 
.0267 



0.137 
.109 
.102 
.091+9 
.0729 
.0325 
-.0292 
-•ilU 
-.123 
-.112 

-.0675 
-.0230 
-.0110 



1.000 

.969 
.872 
.715 
.516 
.288 
.01+63 
-.197 
-.459 

-.81+0 
-.958 

•1.000 



0.126 



.2E 



1.07} 
1.083 
1.113 
1.157 
1.160 

1.178 
1.122 
1.01j2 
.984 
.928 
.906 

.895 



= 1.0928 
To =0.0357 

f „ = 0 

«o 

* 180" 



Pl = 1.0929 

= 0.0263 

IVi = 0 



ISO** 
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01 



TABLE VIII.- INVERSE CMF METHOD FOR 



a = 1 CAMBER LINE VELOCITY DISTHIWTIOI 

NATUNM. mm 



(radians) 




ay,. 


dCx^ 


d5y^ 




&7l 


IT 




'^l 


k 


[cj ='0.67) 


,v 

(Modified 
distribution 
for „ 
Cj =1: 0.8) 


I 

9 

10 


0 


-0.023lf 


-0.0517 


0 


0 


0.0825 


0.0556 


0 


0 


0.914.0 


1.171 


1.194 


-.0122 

-.0194 
-.0210 
-.0161 


-.0192 

-.0101 
.0012 
.0111 


-.0575 
-.0180 
.0069 

0 / 

.0206 


.0286 
•oUoi 

.01;^ 5 
.0209 


.OiliB 
.0287 

.0581 
.0^03 


.0805 
.0710 
.05% 
.0375 


.0585 
.0457 
.0259 
-.Q077 


-.0197 
-.0I169 
-.06k9 
-.0769 


-.2^5 

"kn 
-.672 
-.829 


.820 


1.181 
1.183 
1.185 

I.I03 


1.'205 
1.208 
1.214 
1.222 


11 

12 

11 

15 


-.0075 

-.0012 
.0018 
.01x4 

.0217 


.0154 

,011+1 

.01511 

.0166 
.0076 


.0552 
.0125 
.0206 
.OI1.70 
.0256 


.0031 
-.0056 

.OI3I; 
-.0150 
-.0506 


.0351 
.0127 
-.0152 
-.0525 
-.0569 


.0175 
.0016 
.0037 
.0219 


-.0511 
-.1025 

-.0963 
-.0381 
-.0017 


-.07i;a 

■ Ok79 
.079U 
.077I; 


-.937 
-.992 

-985 
-.902 
-.7li7 


.521 
.103 
.171 
.475 
.711 


.988 

.BlfO 
.843 


1.210 

1.177 
.700 
.798 
.815 


i6 
11 


.0151 

0 


-.0072 

-.0204 
-.0256 


-.Oljl 


-.0580 
-.0382 
0 


-.0317 
-.0183 
0 


.0611I 

.0753 

.0802 


.057U 
.0636 
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